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Abstract 

This paper presents a new mathematical model of vehicular traffic, based on the methods of 
the generalized kinetic theory, in which the space of microscopic states (position and velocity) 
of the vehicles is genuinely discrete. While in the recent literature discrete-velocity kinetic 
models of car traffic have already been successfully proposed, this is, to our knowledge, the first 
attempt to account for all aspects of the physical granularity of car fiow within the formalism 
of the aforesaid mathematical theory. Thanks to a rich but handy structure, the resulting 
model allows one to easily implement and simulate various realistic scenarios giving rise to 
characteristic traffic phenomena of practical interest (e.g., queue formation due to roadworks 
or to a traffic light). Moreover, it is analytically tractable under quite general assumptions, 
whereby fundamental properties of the solutions can be rigorously proved. 

Keywords: traffic granularity, generalized kinetic theory, vehicle interactions, stochastic 
games. 
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1 Introduction 



Modeling vehicular traffic at the kinetic scale offers several advantages in terms of capturing some 
of the relevant aspects of the complexity of car flow. The kinetic approach is indeed suitable for 
an aggregate representation of the distribution of vehicles, not necessarily focused on single cars, 
while still allowing for a detailed characterization of the microscopic vehicle-to-vehicle dynamics, 
which are ultimately responsible for the large-scale behavior of the system. In this respect, kinetic 
models are competitive over microscopic ones, see e.g., [12, 13, 18]. In fact, they require a lower 
number of equations, which makes them more accessible by computational and mathematical 
analysis as far as the study of average collective trends is concerned. On the other hand, unlike 
the prevailing approach at the macroscopic scale, see e.g., [6, 15, 16], kinetic models are not based 
on a priori closures of the equations through fundamental diagrams. In particular, if necessary 
they can recover the latter a posteriori out of a lower-scale modeling of microscopic interactions. 
Actually, it is worth pointing out that not all macroscopic models rely strictly on fundamental 
diagrams. Some models get rid of them by envisaging closures of the balance equations oriented 
to a phenomenological interpretation of the driver behavior [1, 4]. For a more comprehensive 
overview of vehicular traffic models at all scales, the interested reader is referred to [2]. 
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Classically, the kinetic representation of vehicular traffic along a one-way road is provided 
by a distribution function / over the mechanical microscopic state of the vehicles. The latter is 
identified by the scalar position x G and speed v € D^, where Dx,Dy C R are the spatial 
and speed domains, respectively. While the former may cither be a bounded interval, such as 
[0, L], L > being the length of the road, or coincide with the whole real axis, the latter is 
normally of the form Dy = [0, Vmax], where Knax > is the maximum speed allowed along the 
road, or possibly Dy = [0, V^max]j with V^^^ > Knax corresponding to the maximum average speed 
attainable by a single vehicle in free fiow conditions. The distribution function / = f{t, x, v) : 
[0, Tmax] X Dx X Dy ^ R+, where Tmax > is the final time (possibly Tmax = +oo), is then such 
that f{t, X, v)dxdv is the (infinitesimal) number of vehicles that at time t are located between x 
and X + dx and travel with a speed between v and v + dv. 

In the above presentation of the kinetic approach, the spatial position and speed of the vehicles 
are tacitly assiimcd to be continuously distributed over D^ x Dy. However, this does not reflect 
correctly the physical reality of vehicular flow. Indeed, the number of vehicles along a road is 
normally not large enough for the continuity of the distribution function over the microscopic 
states to be an acceptable approximation (like in the classical kinetic theory of gases). Vehicles do 
not span continuously the whole set of admissible speeds; likewise, one cannot expect that they 
are so spread in space that at every point x of the road there may be some. In other words, the 
actual distribution of vehicles in space, as well as that of their speeds, is strongly granular. It is 
reasonable to expect this fact to have a nontrivial impact on the resulting dynamics. It is hence 
worth being taken into account in a mathematical model. 

Recently discrete velocity models have been introduced [3, 7, 10]. The idea is to relax the 
hypothesis that the speed distribution is continuous by introducing in the domain Dy a lattice 
of discrete speeds {vi, V2, ■ ■ ■ , v„}, with vi = 0, ?;„ = V^ax) < Vi+i for alH = 1, . . . , n — 1. 
The physical system is then described by n distribution functions fj = fj{t, x) : [0, Tmax] x 
Dx M+, j ~ 1, . . . , n, such that fj{t, x) = f{t, x, Vj) or, in distributional sense, f{t, x, v) = 
x)Sy.{v). Particularly, fj{t, x)dx is the (infinitesimal) number of vehicles traveling at 
speed Wj, that, at time t, occupy a position comprised between x and x + dx. In this way, the 
granular character of the car flow is at least partially taken into account from the point of view of 
the speed distribution. 

Starting from the discrete velocity kinetic framework developed in the above-cited papers, in 
the present work we undertake the discretization also of the microscopic space variable, in order 
to accomplish the program of a fully-discrete-state kinetic theory of vehicular traffic. The basic 
idea is to partition the spatial domain D^ in disjoint cells (subintervals) such that D^ = U^/i. 
Each cell identifies a portion of the road in which only the number of vehicles is known. No matter 
how rough this assumption may seem, it is nevertheless realistic in view of the spatial granularity 
of traffic, which does not allow one to make finer predictions on the spatial position of vehicles 
within each cell. The microscopic state {x,v) of the vehicles belongs then to the discrete state 
space X {vj}j, which generates a new structure of the kinetic equations not immediately 
deducible from more standard frameworks. Particularly, it is worth anticipating that the resulting 
mathematical structure is not a cellular automaton [8, 11] in spite of the discreteness of the space 
state. In fact, as it will be clear in the next sections, vehicles are not assimilated to fictive particles 
jumping from their current site to a neighboring one with prescribed probability. They actually 
flow through the cells with their true speed, according to a transport term duly implemented in 
the time-evolution equations of their distribution functions. Car-to-car interactions, responsible 
for speed variations, are in turn coded in the aforesaid equations, so that finally the evolution of 
the system is not seen as a stepwise algorithmic update of the lattice of microscopic states. 

In developing the program just outlined, this paper aims at covering the whole path from mod- 
eling to numerical simulations to well-posedness analysis of the resulting mathematical problems. 
In more detail. Section 2 derives the general discrete-state equations, resting on the methods of 
the generalized kinetic theory [2, Sect. 6], to be used as a reference framework for specific models. 
This framework accounts, in particular, for stochastic dynamics of speed variations, according 
to the idea that one-to-one microscopic interactions among vehicles are assimilated to stochastic 
games. Section 3 presents a model obtained from the previous framework, which is then used in 
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Section 4 for producing exploratory numerical simulations of some case studies. These include the 
long-term trend of the system, particularly with reference to the fundamental diagrams predicted 
by the kinetic equations, and the dynamics of queue formation/depletion along a road triggered 
by two different causes: roadworks and traffic light. Finally, Section 5 is devoted to the analyt- 
ical study of some relevant qualitative properties of the initial/boundary- value problems. Under 
quite general assumptions, to be possibly regarded as modeling guidelines, global existence of the 
solutions, uniqueness, continuous dependence on the initial and boundary data, non-negativity, 
boundedness are proved, which provides the modeling framework with the necessary amount of 
mathematical robustness. 



2 Discrete space kinetic framework 

In this section, starting from the discrete velocity setting recalled in the Introduction, we derive 
kinetic structures based on a fully discrete space of the mechanical microscopic states of the 
vehicles. 

As previously anticipated, partitioning the spatial domain in a number of cells of finite 
size is a more realistic way of detecting the positions of the vehicles along the road. Indeed, 
it is consistent with the intrinsic granularity of the flow, which does not allow for a statistical 
description of their spatial distribution more accurate than a certain minimum level of detail. In 
addition, as we will see, it enables one to account easily for some effects due to the finite size of 
the vehicles even in a context in which the actual representation is not focused on each of them. 
To this purpose, a useful partition of the spatial domain is in pairwise disjoint cells U, whose 
union is Dx- For a road Dx = [0, L] of length L > 0, we set 

Dx = [jiu ii,nii,=<d, yii^i2, 

m G N being the total number of cells needed for covering Dx, which depends on the size £i of 
each li . We will henceforth assume that cells have a constant size i, chosen in a such a way that 
L/iGN, thus m = L/£. 

Let fij = fij{t) be the distribution function of vehicles that, at time t, are located in the i-th 
cell with a speed in the j-th class. Then the total number Nij of vehicles in 7j with speed Vj is 
Nij = fiji, while the total number Ni of vehicles in 7j is 

n n 

Notice that pi := Ni/i is instead the local density of vehicles in the i-th cell. By further summing 
over i, one gets the total number N of vehicles along the road: 

m m n 

i=l i=l j=l 

Remark 2.1. If A^max > is the maximum number of vehicles allowed in a single cell, then the 
maximum local density is Pmax := Njaax/^, which defines the road capacity. 

From the fij's, the kinetic distribution function / can be recovered as 

m n 

fit, X, v) = EE/^^W1^.W'^''.(^)' (2) 

1/. being the characteristic function of the cell /j {li^{x) = 1 if a; e /j, lii{x) = if a; ^ /,). In 
practice, / is an atomic distribution with respect to the variable v, like in the discrete velocity 
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framework, and is piecewise constant with respect to the variable x. Particularly, this latter 
characteristic implies that vehicles are thought of as uniformly distributed within each cell. 

Usual macroscopic variables of traffic, such as the vehicle density p, flux q, and average speed 
u are obtained from Eq. (2) as distributional moments of / with respect to v: 

m I n \ m/n \ (4- \ 



Notice that p can also be written as p(t. x) = '^^i Pi{t)li^(x), thus showing that the vehi- 
cle density on the whole road is a piecewise constant function in x built from the local cell 
densities. Analogously, after defining the flux at the i-th cell as qi := J2^=i'^jfijj results 
q{t,x) = qi{t)li,{x). 



2.1 Evolution equations for the /ij's 

A mathematical model of vehicular traffic in a fully discrete kinetic context is obtained by deriving 
suitable evolution equations for the distribution functions fij. To this end, we propose a method 
grounded on classical ideas of conservation laws, which however exploits directly the discrete space 
structure as a distinctive feature of the present setting. 

Fix a cell li and a velocity class j. The basic idea is, as usual, that during a time span At > 
the number Nij of vehicles in that cell with that velocity may vary because of: 

• vehicles leaving the cell /j or new vehicles entering it at speed vj from neighboring cells; 

• interactions among vehicles within the cell 7j inducing speed changes. 

On the basis of these two principles, the following preliminary balance can be written: 

Nij{t + At) - Nij{t) = -N^^\[t, t + At]) + Nlf{[t, t + At]) + Jij{t)iAt + o{At) (4) 

with obvious meaning of the symbols N^^, iVj°"*. The term Jij{t) at the right-hand side accounts 
for instantaneous acceleration/braking dynamics per unit space £, that produce transitions among 
the various speed classes. For the moment we keep it generic, deferring to the next section a 
detailed discussion on its structure. Moreover, the term o{At) indicates that we are neglecting 
higher order effects which may contribute to the variation of Nij during the time At. 

If the time span At is sufficiently small, we can assume that the number of vehicles leaving the 
i-th cell depends solely on the number of vehicles already present in /j at time t. To be definite, 
let VjAt < £ for all j, so that no vehicle can move forward for more than one cell in a single At. 
Then the following form of A^,^"* can be proposed: 

N°;\[t, t + At]) = !^$,,,+iAr,,(i) + o(At), 

where: 

• the ratio VjAt/i accounts for the percent amount of vehicles traveling at speed Vj, which are 
sufficiently close to the boundary between the i-th cell (whence they are coming) and the 
{i + l)-th cell (where they are going) for changing cell during the time span At; 

• the term that we name flux limiter across the interface separating the cells /j, /j+i, 
accounts for the percent amount of vehicles in /j that can actually move to on the basis 

of the free space available in the latter. 

Likewise, the following form of Nl'^ can be proposed: 

N'^{[t, t + At]) = !^$,_i,,Ar,_i,,(i) + o{At), 
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where now the space dynamics involve the (i — l)-th and the i-th cells. 

Plugging these expressions in Eq. (4) and rewriting in terms of the fij 's yields 

fij{t + At)£ - fij{t)£ = i^i,i+,fijit)e - + Jijit)iAt + o(At) (5) 

whence, dividing by £At, taking the limit — 0+, and rearranging the terms produces 

^ + - ^i-i,ifi-ij) = Jii- (6) 

This equation is in principle well defined only in the internal cells of the domain D^, namely 

for 2 < i < m — \. Indeed, at the boundary cells i = 1, z = to the flux limiters <&o.ij ^m,m+i call 
for the two nonexistent external cells Jq, Im+i- In addition, for i = 1 the equation requires also 
the value /oj. In order to overcome these difficulties, we consider that the condition Vj > for all 
j implies a unidirectional (rightward) flux of vehicles, which requires thus a condition on the left 
boundary of the domain. Consequently, the values $o,i and /oj have to be provided as boundary 
conditions. On the contrary, the flux limiter ^m.m+i has to be prescribed at the right boundary in 
order to specify how the natural outflow of the vehicles is possibly modifled by external conditions. 
In standard situations, it will be simply $m,m+i = 1- 

Remark 2.2. The formal derivation of Eq. (6) closely reminds of the standard one of conservation 
laws: the time variation of the conserved quantity is related to some drift plus microscopic con- 
servative interactions typical of kinetic theory. Nevertheless, it is worth pointing out that classical 
procedures, such as the straightforward use of the divergence theorem, for deducing conservation 
laws in differential form cannot be applied in the present context. Indeed, the finite size of the 
space cells is a structural feature of the problem at hand, which prevents one from taking the limit 

Remark 2.3 (Conservativeness of Eq. (6)). The total number N of vehicles along the road, defined 

by Eq. (1), has to balance with the inflow/outflow of vehicles from the boundaries of Dx- We 
can formally check this property out of Eq. (6), after recalling that the conservative interaction 
operator is required to fulflU 

n 

^Jy(t) = 0, Vi = 1,...,TO ,i > 0. (7) 



Multiplying Eq. (6) by £ and summing both sides over i, j gives then 
dN 



= 0, 



.i=l 

whence, summing telescopically in the square brackets, yields the desired result: 

QOi Qm being the incoming and outgoing fluxes of vehicles through the boundaries a; = and x = L, 
respectively. 



2.2 Stochastic dynamics of speed transitions 

The right-hand side of Eq. (6) is the phenomenological core of the kinetic equations, because it 

accounts for the time variation of the conserved quantities fij in terms of microscopic interactions 
among vehicles. The operator Jjj is called interaction operator. It describes modifications of the 
speed Vj of a vehicle traveling in the cell /j due to the action exerted on it by other vehicles. 
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Vehicles arc not purely mechanical particles, for the presence of drivers provides them with 
the ability of taking decisions actively, i.e., without necessarily the influence of external force 
fields. Consequently, the description delivered by Jij cannot bo purely deterministic: personal 
behaviors have to be taken into account, which are suitably modeled from a stochastic point of 
view. Specifically, interactions among vehicles are assimilated to stochastic games among couples 
of players. The game strategy of each player, viz. vehicle, is represented by its speed, while the 
payoff of the game is the new speed class it shifts to after the interaction. Games are stochastic 
because payoffs are known only in probability. In other words, the model considers that vehicles 
might not react always the same when placed in the same conditions due to possible partly 
subjective decisions of the drivers. 

These ideas are brought to a formal Icivel by writing Jij as a balance, in the space of microscopic 
states X of vehicles gaining and losing the state {li, vj) in the unit time: 

if" " \ 

■^y = 2 ^hk{i)Kk{'i')fihfik - fio^V]k{i)fik ■ (8) 

fe=i fe=i / 

This expression is formally obtained by integrating, over a single space cell Jj, the interaction 
operator deduced in previous works dealing with the discrete-velocity kinetic theory of vehicular 
traffic (see e.g., [17]). In this sense, Eq. (8) is consistent with other models relying on similar 
theoretical backgrounds. 
We mention that: 

• Only binary interactions among vehicles are accounted for, which possibly produce a change 
of speed but not of position of the interacting pairs. In particular, it is customary to term 
candidate the vehicle which is likely to modify its current speed vu into the test speed Vj 
(the payoff of the game) after an interaction with a field vehicle with speed Vk ■ 

• The quantitative description of the interactions is provided by the terms r) = {?7hfc(j)}L~fe=i r 
A = {^ft,fc(«)}/,~j'fc=i" „ called the interaction rate and the table of games, respectively. 
The former models the frequency of the interactions among candidate and field vehicles. 
The latter gives instead the probabilities that candidate vehicles get the test speeds after 
interacting with field vehicles (hence Aj^f. is, for j' = 1, . . . , n, the probability distribution of 
the candidate player's payoff conditioned to the strategies Vh, Vk by which the two involved 
players approach the game). Notice that both rjhkii) and ^'^^.(i) may vary from cell to cell. 
Since the table of games is a discrete probability distribution, the following conditions must 
hold: 

0<A^^fe(«)<l, \/h, j, k = l, ...,n, Vz=l, ...,m 

n 

Y^Aii^{i) = l, V/i, fc=l, ...,n, Vi = l, ...,m. (9) 

i=i 

Particularly, from Eq. (9) it follows the fulfillment of Eq. (7). 

• The coefficient £/2 appearing in expression (8) weights the above interactions over the char- 
acteristic length of a single space cell. 

Remark 2.4 (Generalization to nonlocal interactions). The form (8) of the interaction operator 

tacitly implies that candidate and field vehicles stay within the same space cell. However, one can 
also assume that the candidate vehicle interacts with field vehicles located in further cells ahead. 
This amounts to introducing a set of interaction cells 

{It, Ii+i, • • • , Ii+fi{i)} 

extending from the i-th cell up to ii{i) cells in front. The number /x(z) varies in general with the 
index i so as to satisfy the constraint i + fi{i) < m (i.e., interaction cells cannot lie beyond the 
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last cell Im)- Consequently, the structure (8) of modifies as 

^ / i+nii) n i+fJ.{i) n \ 

•^'^■=2 ^ ^ ^h>'^^)^hkC-)fihflkWl- fij X] ^'njk{l)flkWl\ , (10) 
y l=i h,k=l l=i fc=l J 

where are suitable weights accounting for different effects of the interactions based on the 

distance of the Z-th cell from the i-th one. They have to comply with some minimal requirements: 

> y I = i, . . . , i + wi = I Vi = l, ...,m. (11) 

l=i 

Notice that, for = 0, Eq. (11) implies that the only remaining weight Wi is equal to 1, which 
ultimately enables one to recover Eq. (8) as a special case of Eq. (10). 



2.3 Non-dimensionalization 

At this point, it is useful to rewrite Eq. (6) in dimensionless form. To this end, we choose £ 
and I4iax as characteristic values of length and speed, respectively, and we define the following 
nondimensional variables and functions: 

V* — ""^ f* — ff.(t*) — ^ f- ( —^f 

*' max ^ * max \ •'max 

Substituting these expressions into Eqs. (6), (8), we get, after some algebraic manipulations, 

IP* n n 

+ V* (<f>.,,+i/* - <^>i-i,ifu,j) = J2 ^ik{i)Kki^)f:hnk - ftj E '?lfe(*)/4, (12) 

h,k=l fe=l 

where 

£N 

* / .\ -^-^*max /.\ 

VhkW ■■= 757? — VhkW 

^ *^max 

is the dimensionless interaction rate. 

It is worth noting that the new speed lattice is now 

t;* = <■•■<<< <+i <•■•<< = 1, 

while the nondimensional length of the space cells becomes unitary. In addition, one can check 
that the dimensionless local cell density reads 



pKtl := j^P. i^t* = -j;^ L (13) 

-^''max V^^max / -^''max 

whence < p* < 1. Consequently, p* can be directly thought of as a probability density. 

We will henceforth refer to Eq. (12), omitting the asterisks on the dimensionless quantities for 
brevity. 



3 From the general framework to particular models 

Starting from the mathematical structures given in Sect. 2, specific models can be produced by 
detailing the flux limiters, the interaction rate, and the table of games. This amounts to analyzing 
the dynamics of microscopic interactions occurring among vehicles. 



7 



Flux limiters We recall that the term limits the number of vehicles that can actually 

travel across the cells on the basis of the occupancy of the destination cells. 

A prototypical form of such a limiter can be derived via the following argument. The free space 
left in the cell /i+i for vehicles coming from the coll li is Nj^ax—Ni^i. Assuming Ni > TVmax— A^i+i, 
i.e., that the number of vehicles in the i-th cell is greater than such an available space, the percent 
amount of vehicles that can actually move is 



N„ 



— 



H+l 



On the other hand, if < A^max — A^i+i then in the cell /j+i there is, in principle, enough room 

for all vehicles coming from the cell li, thus in this case ^i.i+i = 1- 

Rewriting the flux limiter in terms of the local cell density via Eq. (13), it results 

= { Pi (14) 
if Pi + Pi+i < 1- 

The interpretation is that if the total density in the pair of cells {li, li+i) is greater than the 
(dimensionless) road capacity then the microscopic granularity of traffic starts acting and all 
vehicles cannot freely flow through the cells. Otherwise, the flux limiter has actually no effect. 

Interaction rate The frequency of the interactions can be assumed to depend on the number 

of vehicles in each cell, in such a way that the larger this number the higher the frequency. In 
view of the indistinguishability of the vehicles, a possible simple form implementing this idea is 

Vhkii) = = VoPi, (15) 

'max 

where 770 > is a basic interaction frequency. Notice that this speciflc form does not depend 
explicitly on the speed classes of the interacting pairs. From now on wo will invariably consider 
interaction rates of this kind and we will use accordingly the simplifled notation r]{i). 

Table of games Speed transitions within the space cells depend, in general, on the evolving 

traffic conditions, namely, among others, the local car congestion and the quality of the environ- 
ment (e.g., road or weather). In this paper, we account for the former via the concept of fictitious 
density explained below, which aims at introducing in the table of games the active ability of 
drivers to anticipate the actual evolution of traffic. Conversely, we identify the latter simply by 
a parameter a G [0, 1], a = standing for the worst environmental conditions and a = 1 for the 
best ones. 

The concept of fictitious density, originally introduced in [9], is here revisited in the frame of 
the discrete space structure and of the stochastic game approach to vehicle interactions. We define 
the fictitious car density pi felt by drivers in the cell /j as the following weighted average of the 
actual car densities in the cells li, /j+i: 

f(l for i = 1, m- 1 
Pi-= < f . (16) 

[pm for I = m, 

where ^ G [0, 1] is a parameter related to the anticipation ability of drivers. From Eq. (16) it 

results that if pi^i < pi then pi < pi, thus if the car density in the cell ahead is lower than the 
one in the current cell drivers may be motivated to behave as if the traffic congestion in their 
cell were lower than the actual one. This implies, for instance, a higher inclination to maintain 
the speed, or even to accelerate, in spite of a possibly insufficient local free space, because they 
are anticipating the nearby evolution of traffic. Conversely, if pi+i > pi then it results pi > Pi, 
hence if the car density in the cell ahead is higher than the one in the current cell drivers may be 
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motivated to keep a more precautionary behavior, for instance by decelerating in spite of the real 
local free space. Notice that pi = pi if and only if /3 = or pj+i = pi. In these cases, driver's 
anticipation ability is countered by either a specifically modeled behavior or the instantaneous 
traffic evolution. 

The table of games stemming from these ideas is described in the following. Three cases need 
to be dealt with, corresponding to candidate vehicles h traveling more slowly, or faster, or at the 
same speed as field vehicles k. 

* Vh < Vk- In this case, the candidate vehicle may either decide to maintain its speed Vh 

or be motivated by the faster field vehicle to accelerate; to the speed u/i+i. The better the 
environmental conditions (a), and the more the free room felt by drivers in the i-th cell 
(1 — pi), the higher the probability to accelerate. Such actions are further modulated by 
the actual free room available in the destination cell li+i. This effect is accounted for by 
the flux limiter specifically by the fact that the candidate vehicle can be forced to 

stop (corresponding to the speed transition Vh ^ Vi = 0) as the next cell fills. A technical 
modification is necessary at the boundary ft, = 1 of the speed lattice, for if the candidate 
vehicle is already still then the options of either maintaining the speed or being forced to 
stop by the flux limiter coincide. 



Vh<Vk{ 



Ai,(z) = l-a(l-p,)$i,,+i 
h = l { Al^{i) = a{l - pi)^i,i+^ 
= otherwise 



h>l { 



[Am 

A^,(z) = [l-a(l-p,)]^M+i 
4+i(z) = a(l-p,)^,,,+i 

~ ^ otherwise 



(17) 



• Vh > Vk- In this case, the candidate vehicle may either maintain its speed Vh, if e.g., it can 
overtake the leading field vehicle, or be forced to queue by decelerating to the speed Vk of the 
latter. Specifically, from the expression of the corresponding transition probability reported 
in Eq. (18) below, it can be noticed that this latter event is fostered by bad environmental 
conditions as well as by a reduced free room felt by drivers in the i-th cell. In particular, if the 
candidate vehicle interacts with a stationary field vehicle (fc = 1), the transition probabilities 
are deduced by merging, like before, the first two cases of the expression valid for /c > 1. 



Vh>Vk< 



k > 1 





= 1 - a{l - p^)^i,i+l 


< Am 


= a{l - ft)$»,i+i 




= otherwise 


[Am 


= 1 - 


, Am 


= [1 - a(l - Pi)]^i,i+i 


Am 


= a{l - pi)$i_i+i 


Uikii) 


= otherwise 



(18) 



• Vh = Vk- Finally, in this case the candidate vehicle can either be motivated to accelerate to 

the speed v^^i or induced to decelerate to the speed by the field vehicle; alternatively, 
it may also decide to maintain its speed Vh- The corresponding transition probabilities are 
built from arguments analogous to those discussed in the previous points, thus invoking 
again the concepts of environmental conditions, free room felt by drivers in the i-th cell, and 
modulation effect by the flux limiter. At the boundaries of the speed lattice some technical 
modiflcations are needed, for some output events may coincide while others may not apply 
(e.g., deceleration and acceleration are not possible in the lowest and highest speed classes. 
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respectively). 



Vh = Vk< 



h = 2 



h=l 




Uii« 

^22(0 

^ ^22(*] 
^22(0 
.^22(*] 



= 1 - + (1 - a)/5i$i,i+i 

= [1 - a - (1 - 2a)pi]^i,i+i 



= 1 - a(l - pi)^i,i+i 

= a(l - Pi)^i,i+i 
~ otherwise 



= otherwise 



= 1 - <i>u+i 
) = (1 - a)p^$^,,+i 
= [1 - a - (1 - 2a)pi]^i,i+i 
) = a(l - 



(19) 



2 <h<n 



= otherwise 



h = n 




n-l 
nn 



'^) 



= 1 - 

) = (1 - a)p,$»,,+i 
= [1 - (1 - a)pi]$i,i+i 



= otherwise 



4 Computational analysis 

The vahdity of the model presented in the previous sections, and indirectly also that of the method- 
ological approach which generated it, can be assessed through exploratory numerical simulations 
addressing typical phenomena of vehicular traffic. 

4.1 The spatially homogeneous problem 

A first common case study refers to car flow in uniform space conditions, i.e., when one assumes 
that cars are homogeneously distributed along the road. As a matter of fact, this ideal scenario 
is the tacit assumption behind the empirical study of the so-called fundamental diagrams out of 
experimental data measured for real flows. Such diagrams relate the density of cars to either their 
average speed or their flux, this way providing synthetic insights into the gross phenomenology of 
car flow expected in stationary conditions. The representation at the mesoscopic (kinetic) scale 
is particularly suited to simulate, by mathematical models, this kind of information. In fact, the 
distribution function allows one to obtain genuinely the required statistics (cf. Eq. (3)) from the 
modeled microscopic interactions, rather than relying on them for modeling directly the average 
macroscopic dynamics. Speciflcally, considering the basic assumption of spatial homogeneity, the 
equations to be solved are: 



which are derived from Eq. (12) with fij = fi+ij for i = 1, . . . , to — 1 and all j, along with the 
assumption that the interaction rate does not depend on h, k, cf. Eq. (15). In view of Eq. (14), 
also the flux limitcrs turn out to be independent of i, which makes the transport term at the 
left-hand side of Eq. (12) vanish. Ultimately, the index i of the spatial cell has been removed 
because it is unnecessary in the present context. 

Owing to Eq. (9), Eq. (20) is such that the car density p = fj is conserved in time: 




(20) 
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Figure 1: From left to right, fundamental diagrams of the average flux, the average speed, and the 
speed variance computed from the stationary solutions of the spatially homogeneous model (20) 
for various values of the parameter a indicated in the key. 



therefore one can fix p S [0, 1] at the initial time and then, by means of Eq. (20), look for the 
asymptotic statistical distribution (if any) of microscopic speeds corresponding to the very same 
car density. Following this procedure, the fundamental diagrams of Fig. 1 have been numerically 
obtained for various values of the parameter a representing the environmental conditions in the 
table of games (cf. Eqs. (17)-(19)). It is worth noticing that the other parameters ryo (cf. Eq. (15)) 
and P (cf. Eq. (16)) play instead no role in shaping the asymptotic distributions, thus also the 
fundamental diagrams. Indeed, the first one can be hidden in the time scale whereas the second 
one drops because in spatially homogeneous conditions the fictitious car density coincides with 
the actual one. 

Diagrams in Figs, la, lb have been computed by duly adapting to the present case the formulas 
reported in Eq. (3). They demonstrate that the model is able to catch qualitatively the well-known 
phase transition from free to congested traffic flow [14]. The trend of the average flux q (Fig. la) is 
basically linear for low density, then it becomes markedly nonlinear when p increases, i.e., the road 
progressively clogs up. Notice that for a = 1, i.e., in the ideal situation of optimal environmental 
conditions, the average flux is exactly linear in the free flow phase, which persists until p = 0.5. 
The average speed u (Fig. lb) is initially almost constant and close to the maximum possible 
value u = l (actually, for a = 1 the model predicts that it is exactly constant to 1 until p = 0.5), 
then it drops steeply to zero when the density enters the congested flow range. The free flow 
phase reduces as the environmental conditions worsen. For a decreasing from 1 to 0, it eventually 
disappears completely (cf. the curves for a = 0.5). 

The curves in Fig. Ic show the asymptotic trend of the speed variance 9, defined as 

1 " 

vs. the car density p. They can provide useful quantitative information about the aforesaid phase 

transition, considering that they clearly show that there is a value of p (depending on a) for 
which 9 is maximum. By understanding such a maximum as the signal of a strong change in 
the microscopic behavior of vehicles with respect to the estimated average one, it is possible to 
interpret the corresponding density as the critical density pc separating the free and congested 
regimes. For instance, experimental measurements on the Venice- Mestre highway [5] indicate that 
a realistic range of the critical density is pc S [0, 0.15]. According to the presented simulations, the 
corresponding realistic range of the environmental parameter is a G [0, 0.61]. Of course, different 
roads can produce different empirical evaluations of pc, however the tuning procedure of a based 
on the fundamental diagram of 9 remains conceptually the same. 
A couple of comments are in order: 
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Parameter 



Value 



Description 



m 



n 



a 



10 
6 
1 

0.61 



Number of space cells 
Number of speed classes 

Coefficient of the interaction rate, see Eq. (15) 
Constant environmental conditions 




if j < 5 

%) if 6 < i < 9 
if i = 10 



Cell-to-cell variable environmental conditions 







Coefficient of the fictitious density, see Eq. (16) 



Table 1: Model parameters for the roadworks problem 



Parameter Value Description 



m 10 Number of space cells 

n 6 Number of speed classes 

r]Q 1 Coefficient of the interaction rate, see Eq. (15) 

a 0.55 Environmental conditions 

/3 1 Coefficient of the fictitious density, see Eq. (16) 



Table 2: Model parameters for the traffic light problem 



• the simulations based on Eq. (20) indicate that for a < 0.5 the speed variance 9 attains 
its maximum always at p = 0, i.e., the free flow phase is absent. Therefore the trend of 
the variance does not provide enough information for identifying univocally the value of a 
corresponding to a null critical density, this particular case requiring presumably a more 
accurate investigation; 

• the same simulations indicate also that in no case the critical density is expected to be greater 
than pc = 0.5 (which is the value found for a = 1). This can be explained considering that 
for p > 0.5 cars are, in average, closer to one another than their characteristic size (i.e., the 
average distance between two successive cars is lower than the typical car size) , which might 
prevent any of them from being unperturbed by cars ahead. 

4.2 The spatially inhomogeneous problem 

A second case study concerns the evolution in time of the car density along a road. Nonuniform 
space conditions arc now explicitly considered, the interest being also in the transient trend of 
the traffic flow. A qualitative validation of the model can be performed upon assessing its ability 
to reproduce spatial flow patterns observed in real conditions. Here we focus specifically on the 
formation of queues, presenting two prototypical scenarios which may be of practical interest for 
traffic management. All relevant parameters used in the next simulations arc summarized in 
Tables 1, 2. 

The first scenario refers to queue formation caused by a worsening of the driving conditions, 
for instance a narrowing of the roadway due to roadworks. Using the full Eq. (12), this can be 
simulated by assuming that the parameter a in the table of games varies from cell to cell, a = ai, 
i — I, . . . , m. In particular, inspired by the reasonings proposed for the spatially homogeneous 
problem, we consider a gradual decrease of a from 0.61 at the beginning of the road to 0.5 at its 
end. In the case of the already mentioned Venice- Mcstre highway, this corresponds to a change in 
the expected flow conditions from pc = 0.15 to Pc = 0. Figure 2 clearly shows the formation and 
backward propagation of a queue for such a variable a with respect to the basically unperturbed 
flow of vehicles produced when a is kept constant to its best possible value (a = 0.61) all along 
the road. 

The second scenario refers instead to the alternate formation and depletion of a queue behind a 
traffic light, which, as Fig. 3 demonstrates, are well reproduced by the model. In more detail, the 
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t = 20 



f= 100 



t = 200 



constant a 



variable a 




0.1 0.2 0.3 0.4 0.5 06 0.7 08 09 10 01 0.2 03 0.4 05 06 0.7 08 0.9 



Figure 2: Comparison, at the same computational times, of the traffic flow under constant and 
variable road conditions due to roadworks. The road is initially empty, while a constant density po 
of incoming vehicles, uniformly distributed over the speed classes, is imposed at the left boundary 



at all times. The flux limiters $0,1, 'I'lo,!! at the left and right boundaries are set to 



i-pi 

Po 



and 



to 1, respectively. The dotted line spanning the space cells is the graph of the function a = a{x) 
used in these simulations. 



t = o 



t= 10 



t= 100 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3: Evolution of a queue during several red-green cycles of a traffic light placed at a; = 0.5. 
No cars enter the domain from the left boundary for all t > 0. The initial condition is the pre- 
formed queue illustrated in the left snapshot. The traffic light period is of 20 computational time 
units. 



action of the traffic light on the flow of vehicles is easily simulated, in our framework, by identifying 
the traffic light with one of the flux limiters ($5,6 in the case of Fig. 3), to which a periodic time 
evolution is assigned mimicking the successive red-green cycles of the light. Specifically, during 
the red phase the selected limiter is set invariably to zero, since no car is allowed to fiow beyond 



13 



the traffic light regardless of the actual free space in the cell ahead. Conversely, during the green 
phase the selected limitcr switches to the expression given by Eq. (14), hence the flow of cars is 
again possible according to the usual rules. The initial condition for this problem is, as shown by 
the first snapshot of Fig. 3, a fully formed queue (i.e., p, = 1) standing at a red light with empty 
road ahead. When, at t > 0, the light turns to green, cars in the cell I5 = [0.4, 0.5] can start 
moving only if they are able to see the empty space in front of them, namely if they react to a 
fictitious density p5 lower than the real one = 1. Indeed, by carefully inspecting the table of 
games reported in Eqs. (17)-(19), it results that standing vehicles (i.e., vehicles in the speed class 
v\ = 0) in a full cell (i.e., one with pi = 1) have no chance to shift to the velocity class U2 > if 
the transition probabilities depend on the actiial coll density. Therefore, the concept of fictitious 
density as introduced in Eq. (16) is important for conferring realism to the model. 



5 Qualitative analysis 

The aim of this section is to study the well-posedness of the following class of initial/boundary- 
value problems, generated by the mathematical framework presented in the previous sections: 



dt 

fim = f?j 



which we specifically consider in the mild formulation 



,h,k=l 



(21) 



Jo 



(,s)/,_i,,(s)-$,,,+i[f](s)/i,-(s)) 



Gij[i,i]is) + fijis)Lij[i]{s)}ds 



(22) 



,/o,W-/,(i) 
$0,1 = $(i) 



where 



G.,[f,f]=r7(z)[f] E KkimMk, 

h,k=l 



(23) 



are the so-called gain and loss operator, respectively (cf. Eq. (8)). Notice that unlike the previous 

sections here we have emphasized the dependence of the flux limiters, the interaction rate, and 
the table of games on the set of distribution functions f = {fij}i,j- This notation will be indeed 
useful in the next calculations. 

We organize the qualitative analysis of Problem (21) in two parts: first we prove the uniqueness 
of the solution and its continuous dependence on the initial and boundary data; next we state 
existence. As a preliminary to this, we introduce the functional framework in which we will 
consider Problem (21) and the assumptions needed for achieving the proofs. 

We denote by Xt^^^ = C([0, Tmax]; M™") the Banach space of vector- valued continuous func- 
tions u = u{t) : [O, Tmax] K™", u(i) = (uii(i), . . . , «i„(i), . . . , Umi(i). • • • . Umn{t)), equipped 
with the uniform norm 



sup ||u(t)lii= sup EEl"yWI' 
te[o,Tmax] te[o 
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where || • ||i is the classical 1-norm in M™". We also introduce its subset 

B = |u e : < Uij{t) < 1, J2uij{t) < 1, Vi, j, Vt e [0, T^ax] 

which is readily seen to be closed. We point out that, by virtue of the equivalence of all norms 
in finite dimension, a different p-norm may be considered in M™". However, the 1-norm takes 
advantage more directly of the interpretation in terms of number density: if u G B then ||u(t)||i 
is the density of the distribution u at time t (cf. Eq. (1)). 

A solution to our initial/boundary- value problem is defined as follows: 

Definition 5.1. A function f = f{t) : [0, Tmax] MJ"" is said to be a mild solution to Prob- 
lem (21) if f e S and f satisfies (22). 

Notice that i € B guarantees that f fulfills the basic requirements of a physically acceptable 
distribution function: non- negativity and boundcdncss in each spacc/spccd class in the whole time 
interval of existence. 

The following box summarizes the hypotheses, to be possibly regarded as modeling guidelines, 
under which we will achieve the proofs of the announced results. 



Assumptions for Problem (21) 

1. The initial and boundary data f^j, fj belong to B. 

2. The fiux limiter [f ] is defined by a function $ : [0, 1] x [0, 1] [0, 1] acting on 

the densities of the two adjacent cells Ji, /j+i such that $j^i_|_i[f] = $(pi, Pi+i)- The 
function $ has the following properties: 

(i) < $(m, v) < 1, Vu, V e [0, 1]; 

(ii) $(«, v)u < 1 - w, V(w, w) G [0, 1] X [0, 1] such that u + v>l; 

(iii) |$(U2, V2) - ^(wi, vi)\ < Lip($)(|'U2 - Ui\ + \v2 - wi|), V(mi, vi), iu2,V2) e [0, l]x 
[0, 1]. 

3. The boundary datum $ is as in Assumption 2, i.e., $ = $(p, pi) where p = J2]j=i fj 
the prescribed density of incoming cars. 

4. The interaction rate r]{i) is: 

(a) non-negative and uniformly bounded, i.e., there exists a constant ^ > such that 

< 7?(i)[u] < Vz = 1, . . . m, Vu G S; 

(b) Lipschitz continuous with respect to the distribution in B, i.e., 

\v{i}[um - v{i}[vm\ < Lip(r?(i))||u(i) - v(t)||i, G [0, T^^], Vu, v G B. 

5. The elements A-f^i.{i) of the table of games satisfy Eq. (9) and the following Lipschitz 
continuity condition: 

AL(i)Mt) - Ai,(i)[vm < Lip(4,(i))||u(t) - v(i)||i, G [0, T^ax], Vu, v G B. 



Wc henceforth assume that Assumptions 1 4a arc always satisfied, indeed they are needed in 
all of the next results. On the contrary, we will explicitly mention Assumptions 4b, 5 only when 
they will be actually used. 
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5.1 Uniqueness and continuous dependence on the data 

The following theorem states that the mild solution to Problem (21) is unique in B and that it 
depends continuously on the initial and boundary data. 

Theorem 5.2 (Uniqueness and continuous dependence) . Under Assumptions 4b, 5, let{i'^, f, 3>-^), 
(g°, g, §^) be two sets of initial and boundary data and f , g e S two corresponding mild solutions 
to Problem (21). Then there exists C > such that: 



l|f-g||oo<C 



llfO „0 



g 111 



\^f{t)-^^{t)\ + \\m-m\i) dt 



In particular, there is at most one solution corresponding to a given set of initial and boundary 
data. 

Proof. Subtracting term by term the mild equations satisfied by f, g (cf. Eq. (22)) and then 
summing over i = 1, . . . , m and j = 1, . . . , n gives: 



iifW-gWiii<iif°-g°iii+EE 

i=i j=i 



/ kj(^i-i,i[f](s)/i-i,i(s) -^i-i,i[g](s)5i-ij(s) 
Jo 

+ $i,i+i[g](s)5y(s) - $i,i+i[f](s)/ij(s))| ds 
(|G.,[g, g](s) - G„-[f, f](s)| + \fijis)Lij[f]is) 

- 9ij{s)Lij[g\{s)\^ ds . 



By exploiting the Lipschitz continuity of the various terms and extracting from the sums the 
boundary terms corresponding to i = 1 we then discover: 

rn n 

EE l*»-l.iW(s)/i-l,j(s) - $i-l,i[g](s)5i-lj(s) + $M+l[g](s)5ij(s) - ^i,i+l[f]{s)fij{s)\ 

i=l j=l 

< 2(2Lip($) + l)||f(s) - g(s)||i + ir (s) - $5(s)| + ||f(s) - g(s)||i, (24) 



J2J2\(^iAs,s]{s)-Gi^[f,f]{s) 



--1 = 1 



< 



+ E I Lip(?7(z)) + r? E Lip(A^^fe(^)) 

h,k,3=l 



f(s)-g(s)||i, (25) 
(26) 



Y.Y.\^iAs)^iAm - 9i3{s)LiAE\{s)\ < 2iyn + n2ELip(»?W) ||f (s) - g(s)||i, 

i=l j=l \ i=l / 

whence 

\\m - g(t)||i < ||fO - gO||i +C f (||f(a) - g(5)||i + \¥{s) - $«(s)| + ||f(a) - g(a)||i) ds 

Jo 

where C is a suitable constant deduced from all those appearing in Eqs. (24)-(26). Gronwall's 
inequality implies now: 



m-m\i<e 



ct 



„0| 



llf - 



+ / (|$/(s)-$'^(s)| + ||f(s)-g(s)||i) ds 
Jo 



whence, taking the suprcmum over t £ [0, Tmax] on both sides, we obtain the continuous depen- 
dence estimate. Uniqueness follows straightforwardly from it by letting f" = g°, f = g, and 

$/ = $9. □ 
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5.2 Existence 



In this section wc prove the existence of the mild solution to Problem (21). The idea of the proof 
is to consider at first the model at discrete time instants = nAt^, where the index n labels the 
discrete time and the index k is a mesh parameter denoting the level of refinement of the time grid. 
In particular, the time step At^ is chosen such that it tends to zero for k oo. After establishing 
some technical properties of the iterates of such a discrete-in-time model (cf. Lemma 5.3), we 
pass from discrete to continuous time by interpolating the iterates and we prove that in the limit 
k ^ oo this interpolation converges to a distribution function f € B (cf. Lemma 5.4). Finally, we 
show that f is precisely a mild solution of to Problem (21) (cf. Theorem 5.5). 



Step 1 We consider Problem (21) at discrete time instants, cf. Eq (5): 

+A4(G.,[f"•^^'^]-/:)%■[f">^]). (27) 

The following lemma gives some properties of the iterates /"j^- 
Lemma 5.3. If At^ is sufficiently small then f"''* € B for all n, k > 0. 

Proof. The proof is by induction on n, assuming that f"''* € B and that the initial and boundary 
data have been chosen consistently with the properties of the distribution functions in B. 

(i) N on- negativity of the iterates. Taking into account the non- negativity of the terms containing 
$i_i,i[f"''^] and G^J[f"'^ f"^''], it results: 

f:;''' > f^j'[l - At,(^,$,,+i[r.^] + L,,[r.^])] > f^j'[L - At,{l + L,,[f".^])], 
being vj, ^i^i+ilf'^] < 1. Moreover iy[f"'''] < fjp"' < fj, hence if Atk < 1/(1 + ??) we deduce 

f:;'''>f::^'[i-Atk{i+fj)]>o. 

(ii) Boundedness of the iterates. Considering that Vj < 1 wc have: 

f",;'^' < f^t + Atk($i-i,i[r'1/^i,, + G,,[r'\ f"'^] - /^'%,[r'1). 

But 



thus we conclude 



n,ki j-rijk 



(28) 



Now, if f"-^ = 1 then from Assumption 2(ii) it follows f"^^'^ < 1 as desired. Conversely, if 



f"-^ < 1 then according to Assumption 2(i)-(ii) it results 



1 II 



n.k 



• r n.k I n,k ^ 
n.k ^ i n.k -r n.k . n.k ^ 

PiU < i-Pi' if Pi' +PiU < 1 



Hence we can continue the estimate (28) as: 



<l-pf. 

(29) 



f^;'''<f^j'+At,\i-,^. 



n,k 



</5' + A(ai + 2f;)(l-/5") 



n,k\ 



whence we see again that f.+^'^ < 1 if Aik < 1/(1 + 2fj). 



17 



(iii) Boundedness of the sum of the iterates. In order to prove that Y^^=i fit — ^® note that 

from Eq. (9) it follows Ej=i(Gy[f"'\ f"'''] - fij^UA^"'^]) = 0. Consequently, by Eq. (29) 
we obtain the thesis, since 

E fT'^ < pf + Ai,$,_i,[r'V."i\ < pf + Aik(i - pf) < 1 

provided Af^ < 1. 

Finally, we note that all properties proved so far hold simultaneously if Atk < 1/(1 + 2?7). □ 

Step 2 Now we pass by interpolation from discrete to continuous time. To this end, we introduce 
the following function f = P(t) : [0, Tmax] — >■ K"*" interpolating piecewise linearly the iterates 

J ij ■ 



n = l 



At. r^' At. 



Ite (30) 



Here A^k is the total number of time steps in the k-th mesh. We assume that A^k and Aty are 
chosen in such a way that A^kAtk = Tmax independently of the refinement parameter k. 
For the function P we can prove the following result. 

Lemma 5.4. There exists f G B such that 

lim ||f^-f||oo = 0. 

k— >cx) 



Proof. We first claim that P € S. In fact, for all t £ [0, Tmax] there exists n < A^k such that 
t e [iLi, t^l thus: 



V Atk ' Atk 

and according to Lemma 5.3 this value satisfies all constraints required pointwise in time to the 
components of the functions in B. Thus the claim follows from the arbitrariness of t. 

Next, our goal is to show that the sequence {P}k>o is relatively compact in Xt-^^^- Indeed in 
such a case we can conclude that (up to subsequences) it converges uniformly to some f G ^T„ax • 
Moreover, since {P}k>o C B and B is closed, we also have that f G S. We prove the relative 
compactness by means of the Ascoli-Arzela compactness criterion: 

(i) {f'*}k>o is uniformly bounded. Indeed, since P € S we can write 



m n 

< m < 00. 



sup 1 1 f 1 1 oo = sup sup V V fij {t) 
k>0 k>OtG[0,T„J^p^ I 

(ii) {P}k>o is equicontinuous. In fact, let ti € t%] and t2 G [t^-i, t^] with m < n (i.e., 

ti <t2), then we can write: 



n-2 



= [i^(i2) - fUil-i)\ + [f^it) - i^(ii)] + E(/^''' - 4') 

=: A1+A2+ A3. (31) 
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We begin by estimating the term ^43. By Eq. (27) and applying the same arguments as in 
the proof of Lemma 5.3 we have: 



< 2(1 + f])Atw, 



thus 



fi-2 



l^sl < E W'' - fit < 2(n - 1 - m)(l + fj)At^ = 2(1 + fj){tU " 4)- (32) 



As for the term A2, we observe that taking into account Eqs. (30), (32) we can write 



\Ao\ = 



fijit^) - fij{tn-i) = fij{t2)-fij 



n-l,k 



Likewise, 

i^ii= fij{ti)-fUt^) = f^'^-kit,) 



/»n,k /-n— l,k 
J ij J ij 



<2(l + i?)(f2-fti)- (33) 



1 



ti-t\ 



m-l 



/•m,k /•m — l,k 
hj ~ J ij 



<2(l + »?)(4-ii)- (34) 



Finally, by plugging Eqs. (32) (34) into Eq. (31) we get 

fij{t2) - fiji^ < 2(1 + ry)(4 -h+h- tti + iti - 4) = 2(1 + f]){t2 - h), 
whence 

\\tit2)-i^{ti)\\i<2mn{l + fj){h-ti). 

This implies the desired equicontinuity of the family {f'^}k>o- In fact, for every e > it is 
sufficient to take 1^2 — ill < 2mn(i+fj) ' independently of k, in order to have ||P(t2)—f'*(ii)||i < 

Step 3 Now wc arc in a position to state that a mild solution to Problem (21) does indeed exist. 

Theorem 5.5 (Existence). Let Assumptions 4b, 5 hold. Then f G B found in Lemma 5.4 is a 
mild solution to Problem (21) in the sense of Definition 5.1. 

Proof. Owing to the convergence result stated in Lemma 5.4, we can guess that, for a fixed k, 

the function f'^ defined componentwise by Eq. (30) is a sort of approximation of a possible mild 
solution to Problem (21). Therefore, if we plug /]^- into Eq. (22) we expect a reminder e^^ to 
appear: 



fij(t) -fii+l {«.(*M+i[fl(^)/];(«) - ^i-iAi'Ks)fUj 

-G,,[f^ i%s) + fl{s)L,,[f']{s)} ds = j\%{s)ds. 



(35) 



At this point we would like to take the limit k 00 using Lemma 5.4. This will tell us that 
the distribution function f found in the latter is indeed a mild solution to Problem (21) provided: 

(i) the left-hand side of Eq. (35) tends to the corresponding expression evaluated for f rather 
than for P; 
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(ii) the right-hand side of Eq. (35) tends to zero. 

In order to prove (i) it is sufficient to observe that all terms appearing in the integral are 
bounded by an integrablc constant, thus by dominated convergence it is possible to commute the 
limit in k with the integral in t. Using Lemma 5.4 and taking into account Eqs. (24)-(26) with 
g = P yields then: 



lim 

k— foo 



Mt) - /° + /* {V^{<^i,i+l[f]{s)fi,{s) - $,_i,,[f](s)/._i,,(s)) 
Jo 

-Gij[f, f]{s) + Ms)Lij[f]is)} ds 

as desired. 

It remains to prove (ii), for which we need to find an expression of the term ej-^-. Using Eqs. (27), 
(30) we write: 

1 fk Nk 
n=l 

+Gij[r-''\ f"-i'^] - /;}-^'%,[f"-i'^]} ; 



M=in=i I mk V / 

I ^ ~ ^n-l rn-l,k ( rn,V. _ j^n-l,k\ 

^ Atk V''' ) 



n Wk 



■f k r rpkl FfklJ rn-l,krn-l,k . ^ ^n-1 j^n-l,k / j-n,k rn-l 

4-^ii[f ] = Z^Z^'^ff ]S 4- fik + fij [fik - fik 



fc=l n=l 



+ ■ 



t-t 



Aik 



Aik 

n — 1 ^n — l.k / j'n,k /-n — l,k\ 
r J ik yJ ij ~ hj J 



then we take the time derivative of both sides of Eq. (35) and insert these expressions at the left- 
hand side. As a result, we get an expression of e^j in which terms featuring the difference f"j — 
fij~^'^ appear along with others involving differences between pairs of flux limiters, interaction 
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rates, and tables of games evaluated at f" ^'^ and P, respectively. Applying estimate (32) to the 
first ones and invoking the Lipschitz continuity of the second ones we finally obtain: 

for a suitable constant C > independent of k, hence 

|e^j.(s)| < CAikT„,ax 

which completes the proof. □ 

CoroUciry 5.6 (Improved regularity). The mild solution to Problem (21) is of class in 
[0, T^ax], hence it is actually a classical solution. 

Proof. The continuity of f implies that the right-hand side of Eq. (21) is continuous in t for each 
i = 1, . . . , m and j = 1, . . . , n. Hence is continuous. □ 
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